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Abstract 

We use the method of dynamical mean field thoery, to study the effect 
of Kondo resonance on optical third harmonic generation (THG) spectra 
of strongly correlated systems across the metal-insulator transition. We 
find that THG signals are proportional to the quasiparticle weight z of 
the Kondo peak, and are precursors of Mott-Hubbard gap formation. 

PACS 42.65.-k; 71.27.+a; 71.30.+h 
Keywords Kondo resonance; Third harmonic genration; Metal-insulator transtion 

1 Introduction 

One of the earliest and most important triumphs of the Dynamical mean field 
theory (DMFT) was to elucidate the nature of metal-to-insulator transition 
(MIT) driven by strong electron correlation. It was demonstrated that this 
transtion is accompanied by a Kondo resonance at the Fermi surface [2] , width 
of which is proportional to the quasiparticle weight z, and approaches to zero 
as the Mott insulating side (U — > {7~) is approached from the metallic side. 
Therefore in nonlinear optical spectra, we expect some signals that vanish in 
the limit U — > , as a result of qusiparticle weight z ^ 0. 

Equivalently one can employ the pressure to tune U c and hence switch from 
Mott insulating side to a regime in which such signals arise due to Kondo peak. 
This is a unique feature of correlated insulators, that has no analogue in band- 
insulators. In case of band insulators, to generate quasiparticle states in the 
mid-gap, one may need to dope them, while in Mott insulators such mid-gap 
states can be created by applying pressure to a Mott insulator which is in the 
onset of MIT. Possible examples would be a class of materials known as Kondo 
insulators 3 . (A typical example of this family is YbBi2 0]). 

2 Method of calculation 

One of the applications of nonlinear opitcal materials is the third harmonic 
generation (THG) [5], that requires materials capable of providing strong enough 
nonlinear optical signals. To study this process in correlated electronic systems, 
we have developed the theory of nonlinear optical response in the context of 
DMFT 6 which is also a suitable theoretical tool to study the effect of Kondo 
resonance in nonlinear spectra. In d — ► oo limit where DMFT becomes exact, 
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vertex corrections to current operators identically vanish [Sj so that the THG 
sucseptibility within DMFT approximation becomes [H] 

X^^—j^^j^j^j^, (1) 

where £j = u + jv — ~Er(uj + jv) + i|Ej(w + jv)\ for j — 0, 1, 2, 3, and D(e) = 
2a/1 — s 2 /tt is the semicircular DOS corresponding to the normalized hopping 
t= 1/2. 

In the above formula, (i) the integration over e corresponds to summation 
over intermediate states in conventional expressions |Hj which are usually used 
for systems with discrete energy levels, and (ii) the matrix element effects are 
encoded in E(a>), real and imaginary part of which have been denoted by Ej? 
and £j, respectively. It is very crucial to note that we have used absolute value 
of the imaginary part of the self-energy. This is indeed a necessary step to 
transform from time-ordered four-current, to fully retarded one0. 

We obtain self-energy £ by solving the iterated perturbation theory equa- 
tions of DMFT for the semi-circular density of states 0. In DMFT approxima- 
tion, this quantity is purely local. 

The critical value of Hubbard U corresponding to semicircular Bethe lattice 
is U c ~ 3.35, above which we are in Mott insulating phase (Fig. [JJ. As can be 
seen in Fig. ^ the sharp Kondo resonance state at Fermi surface in the onset of 
transition to Mott insulating phase, is characteristic of correlation driven MIT. 



3 Results 

Fig. |21 shows the imaginary part of the THG susceptibility x THG ( l/ ) °f e< t- O 
as a function of incident photon energy v for various values of Hubbard U close 
to MIT. Dotted line corresponds to U = 3.5 which is in the Mott insulating 
phase. The Mott-Hubbard gap is manifested by a region in which Imx THG is 
identically zero. The onset of three-photon absorption as demonstrated earlier 

corresponds to the frequency where Imx THG takes up, and is 1/3 of the 
Mott-Hubbard gap (E g ). 

Three-photon resonances in this formulation correspond to the peak A in 
Imx THG around v ~ 0.8 (dotted line), while two-photon features appear as a 
dip B around v ~ 1.7 0- As can be seen in Fig. [5J except for slight shift in the 
location of the dip B, the two-photon features are not much affected by moving 
to lower values of U in U < U c region. In contrast, the three-photon features is 
entirely distorted as we cross the U c and enter the Kondo regime of U < U c . 

Below the critical value U c ~ 3.35 the system is in metallic side, with a sharp 
Kondo peak present at the Fermi surface |2J. For U = 3.3 (dashed line) just 
below the critical value, there appear peaky structures below v ~ 1. Among 
these structures, weak features denoted by A' and B' in Fig. [21 start to appear 
in the THG 'gap' region (0 < v < E g /3). These features seem to correspond to 
three-photon and two-photon rosonances between the lower Hubbard band and 
the Kondo resonance state (schematically depicted in Fig. [JJ, as they vanish 
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Figure 1: Density of states within DMFT at the onset of transition to insulating 
state (dashed line) and in the Mott insulating state (dotted line). A and A' 
represent schematically the three-photon resonance processes in ^ THG of Fig. |3 
Such three-photon resonances appear as peaks A and A' in Fig. [21 as explained 
in Ref. Similarly B and B' in this figure schematically represent two-photon 
resonances denoted by B and B' in Fig. [21 



in z — > limit. Moving further away from MIT, quasiparticle weight z grows 
and hence A' and B' features also start to grow proportional to z. To see this 
proportionality, note that x™ G m ec l- Q) lii dominated by D(e) effect, and 
hence in the limit z — ► one can replace the Kondo peak with a sharp peak of 
strength z. Then decomposing to partial fractions, each with a single — e) 
denominator in eq. 1(1} . one can see that j— photon features for j — 1,2,3 will 
be proportional to z. 

Implication of this observation for Kondo insulators is that, an external pres- 
sures on the order of a fraction of GPa can lead to instability in the functioning 
of optical devices built on them by filling in the gap and hence reducing the 
contrast of the output signal. 
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Figure 2: Imaginary part of X THG across the Mott metal-to-insulator transition. 
Critical value of U is ~ 3.35. Solid line corresponds to U = 3.2, dashed line 
corresponds to U = 3.3, and dotted line corresponds to U = 3.5, which is in the 
Mott insulating side of MIT transition. The onset of three-photon absorption 
for U — 3.5 is where imaginary part of x THG starts to become non-zero (around 
v ~ 0.5). 
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